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ABSTRACT 

Aims. We searched for the presence of extended emission-line regions (EELRs) around low-redshift QSOs. 

Methods. We observed a sample of 20 mainly radio-quiet low-redshift quasars (z < 0.3) by means of integral field spectroscopy. 

After decomposing the extended and nuclear emission components, we constructed [O iii] A5007 narrow-band images of the EELR to 

measure the total flux. From the same data we obtained high S/N (>50) nuclear spectra to measure properties such as [O iii]/H/3 flux 

ratios, Fe ii equivalent widths and Hj8 line widths. 

Results. A significant fraction of the quasars (8/20) show a luminous EELR, with detected linear sizes of several kpc. Whether or not 

a QSO has a luminous EELR is strongly related with nuclear properties, in the sense that an EELR was detected in objects with low 

Pen equivalent width and large Hj8 FWHM. The EELRs were detected preferentially in QSOs with larger black hole masses. There is 

no discernible relation, however, between EELR detection and QSO luminosity and Eddington ratio. 
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1. Introduction 

Quasars may show extended emission-line regions (EELRs) 
stretching over characteristic scales of ten to hundred kpc. 
Luminous EELRs have been found in p articu lar arourid 
lectrum radio-loud qu asars (RLQs, iBoro son & Ob 
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Heckmanetal. 1991; Crawford & Vanderriesti i200| 
iStockton et al. 2002; Fu & Stockton 2006 ) and also radio galax- 
ies (e.g. IVillar-Martin et al.l l2006l l2007h . There are also indi- 
cations that radio-qtiiet quasars (RQQs) can have EELRs (e.g. 



Boroson et alll 19851; [Stockton & MacKentvlll987l;lBennert et atl 



2002). Most existing studies, however, are based either on long- 
slit spectroscopy or on narrow-band imaging; while the former 
technique will generally capture only part of any extended emis- 
sion, the latter does not provide any spectroscopic information. 

Integral field spectroscopy is a relatively new and power- 
ful tool to study AGN host galaxies. The combined diagnos- 
tic power of imaging and spectroscopy allows for a flexible 
treatmeant of the data, facilitating the construction of broad- or 
narrow-band images, of 2-dimensional velocity maps as traced 
by emission or absorption lines, or the spectral analysis of se- 
lected regions within the field of view. 

Here we present preliminary results from a study of 20 low- 
redshift QSOs, most of which are RQQs, observed with an in- 
tegral field unit (IFU). After decomposing nuclear and extended 
emission we identified 8 quasars with clearly detectable EELRs. 
We focus on exploring the relation between the nuclear spectral 
properties of a QSO and the existence or non-existence of a lu- 
minous EELR. A full account of our IFU study will be given 
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in a separate paper (Husemann et al., in preparation; hereafter 
Paper II). 



2. Observations and Data Reduction 

We targeted a sample of 20 low-redshift (z < 0.3) QSOs, se- 



lected from the Palomar-Green Survey (PG, ISchmidt & Green 
'1983) and from the Hamburg/ESO Survey (HES; IWisotzkiet al. 
(2000); Wisotzki et al., in prep.). The QSOs have apparent mag- 
nitudes in the range of V ~ 14 - 17; corresponding to typical 
luminosities of My < -23. Three objects are steep-spectrum 
RLQs, another one is classified as a flat-spectrum RLQ; all the 
other sources are known to be radio-quiet. To our knowledge, 
this is the largest sample of bright quasars observed with an IFU 
so far 

We employed the Potsdam M ulti-Aperture Spectro- 
photometer (PMAS; lRothetal.ll2005h . mounted on the 3.5 m 
telescope at Calar Alto Observatory, for two observing runs in 
September 2002 and in May 2003. The PMAS lens array was 
used at a spatial pixel ('spaxel') scale of 0.5" x 0.5", which 
for 16 X 16 spaxels resulted in a field of view of 8" x 8". At 
the typical redshifts of our targets we were able to sample the 
quasar environment with a resolution of ~ 1 kpc, and cover 
linear extents of up to ~ 20 kpc in diameter 

All observations used the V300 grating to map the spectral 
features around the Ha and H/? emission lines in a single setup. 
This yielded a moderate spectral resolution of approximately 
6.5 A FWHM as measured by night sky emission lines. The to- 
tal integration time for each object was 1-2 hours divided into 
several 30-45 min. exposures. Each individual exposure reaches 
a S/N level for the integrated nuclear quasar spectrum of already 
more than 50 per pixel in the continuum. 

The raw data were homo geneous l y redu ced using the PMAS 
reduction package P3D by iBeckeiJ (|2002|) . Due to substantial 
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fibre-to-fibre variations in the wavelength calibration of more 
than 0.7 A rms, we used elements of the general reduction pack- 
age R3D for fibre-fed IFUs (Sanchez 2006) to reduce this scatter 
to 0.25 A rms. We extracted and subsequently subtracted a mean 
sky spectrum from those spaxels in each science datacube that 
were deemed free of emission from the target itself. The data 
reduction will be described in full detail in Paper II. 

3. Results 

3.1. Quasar-EELR decomposition 

In order to measure the properties of EELRs, any possible con- 
tamination by the quasar nuclear emission had to be minimised. 
Three components are expected that need to be carefully de- 
blended: Nuclear emission from the AGN point source; the stel- 
lar component of the host galaxy; and any possible resolved ex- 
tended line emission. We found that in our luminous quasars, no 
significant evidence for a resolved stellar continuum above the 
noise could be found in our data, and we only needed to deblend 
the nuclear source and the extended line emission. 

After running se veral tests we ado pted the iterative decom- 
position method by IChristensen et al.l (12006.) . Briefly, a quasar 
spectrum is extracted from the spaxels at the position of the 
quasar, assuming that the contribution of the EELR and the stel- 
lar component is negligible for these spaxels. By dividing the 
spectrum of each spaxel by the quasar spectrum in the wave- 
length range bracketing a broad emission line (e.g. Ha or HyS), 
an appropriate scale factor is determined for each spaxel. A com- 
plete datacube of the quasar emission using the previously de- 
termined scale factors is then constructed. By subtracting the 
quasar datacube from the original one, a residual datacube con- 
taining any possible EELR is obtained. The process can be it- 
erated by extracting a new quasar spectrum from the previous 
reconstracted quasar datacube; however, our solution converged 
already after the first iteration. This method showed the best re- 
sults and minimal PSF subtraction artefacts compared to other 
stra tegies such as off-band image subtractio n or PSF synthe- 
sis jJahnke et al.1 120041; ISanchez et"aLl |2004 that we also ex- 
plored. However, the method is not perfect, and there are fre- 
quently some oversubtraction signatures at the quasar position 
(cf. Fig. [Til. This happens because the extracted quasar spectrum 
contains a small fraction of the extended narrow emission line 
even if the scale factors were properly determined by the broad 
lij3 line. 

For the present paper we focus on the narrow wavelength 
range containing the H/? and the [Oiii] doublet /l/i4959,5007. 
The scale factor can therefore safely assumed to be a constant for 
each spectrum. Collapsing the residual data cube within +20 A 
centred on the nuclear [Om] peak wavelength, we could then 
produce narrow band [O iii] images for each quasar. Figure [1] 
presents the results, displaying all 8 quasars with a detected 
EELR in [Om]. In the bottom row we also show the residuals 
for a representative subset of 3 (out of 12) undetected sources, 
to give the reader a visual impression of the significance of the 
detections. 



3.2. EELR fluxes 

Examining the extracted [O iii] narrow band images after sub- 
traction of the quasar nuclei, we found that a significant fraction 
(8/20) of our objects contain a luminous EELR. We measured 
the total fluxes in these extended emission components by sum- 
mation of all spaxels in the narrow band images, except those 
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Fig. 1. Nucleus- subtracted [Om] narrow-band images extracted 
from the datacubes. All 8 quasars with clearly detected EELR 
are shown in the top three rows. The bottom row presents three 
example cases with non-detected EELR for comparison. In each 
panel, the position of the quasar is indicated by a green cross. 
Note that the apparent holes at most quasar locations are due 
to oversubtraction and thus not real (see text for details). The 
colour coded line flux is in units of 10"'^ erg/s/cm^/arcsec^. All 
frames are 8" on the side and oriented such that N is right and E 
is up. The white marker shows a 5kpc linear scale at the quasar 
redshift. 



used to extract the quasar spectrum as these spaxels are slightly 
affected by oversubtraction of the quasar emission. In the case 
of PG 0157+001, the EELR is evidently larger than the PMAS 
field of view (see Fig. [T]!, and a significant fraction of the total 
flux may be missed by our data. We therefore treat the measured 
EELR flux for this object as a lower limit. In all other cases 
we find the EELRs to be reasonably well contained within the 
recorded images. 

For the objects showing no sign of an EELR we could only 
estimate upper limits to fluxes. We simply adopted the extended 
flux from the faintest EELR measured in our sample as the upper 
limit for all those objects where there was no EELR detected. 
This is clearly a very conservative assumption which we plan 
to improve in the future by running dedicated simulations. The 
adopted limit was then scaled to the continuum flux density at 
5 100 A of each quasar continuum to account for the contrast be- 
tween nuclear and extended emission. 

Three of o ur QSOs were already i maged through narrow- 
band filters by IStockton & MacKentyl d 19871) . Despite the dif- 
ferent observing techniques, our [O m] EELR luminosities agree 
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Fig. 2. Left panel: Nuclear [Oiii]/H;6 flux ratio plotted against 
the nuclear Fe ii equivalent width. QSOs with a detected EELR 
are denoted by the open squares, the undetected ones are shown 
as upper limits. The dashed line indicates a suggestive separa- 
tion between the detected and non-detected EELRs. Right panel: 
Total EELR flux, normalised to the nuclear H/S flux, plotted 
against the FWHM of the broad H/3 component. Note the lower 
Umit for PG 0157-1-001 for which the EELR is larger than our 
field of view. 



within a factor of ~ 2-3 with their measurements. Three QSOs 
were furthermore co ntained in the narrow-band HST study by 
iBennert et al.l (|2002|) . which however revealed mainly high- 
surface brightness features less than 1" to the QSO nuclei. 

In order to obtain a measure of EELR strength indepen- 
dent of the spectro-photometric calibration. Galactic foreground 
extinction and scaling with quasar luminosity, we divided the 
EELR fluxes (and upper limits) by the fluxes in the broad H/3 
components of the individual quasars. 



3.3. Spectral properties of the quasars 

As mentioned above, each datacube provided us also with a high 
S/N spectrum of the central quasar over essentially the full opti- 
cal wavelength range, albeit at moderate spectral resolution. We 
extracted the nuclear emission line properties by fitting a super- 
position of simple analytical functions to the line profiles, using 
the Levenberg-Marquardt x^ minimisation algorithm. For our 
narrow wavelength range, a straight line was sufficient to model 
the local continuum of each spectrum. The broad Balmer lines 
were modelled with two or three independent Gaussian compo- 
nents, in accordance with previous experiences by several others 
(e.g. lSulentic et alJl2002HVeron-Cettv et al . 2006). We fitted the 
[O m] AA4959, 5007 doublet as a system of two Gaussians with 
the same line width, fixed rest-frame wavelength separation and 
a fixed line ratio of 1:3. 

Both H/3 and the [Oiii] doublet can be heavily affected by 
Fe II emission line complexes. The most prominent is the Fe ii 
AA4924, 5018 doublet, which we modelled similarly to the [O iii] 
doublet, but with a fixed fine ratio of 1.28 (/i5018/.i4924). We 
also used the Fe ii complex (5 100A-5405A) longwards of [O iii] 
to measure the strength of the Fe ii emission in those cases where 
the Fe ii emission is weak. A linear quasi-continuum was esti- 
mated from the bracketing emission-line free regions and sub- 
tracted from the quasar spectrum. We integrated the residual flux 
in the given wavelength interval to estimate the Feii flux and 
equivalent widths. 



3.4. Black hole masses and accretion rates 

We estimated black hole masses (Mbh) for all objects from 
our single epoch quasar spectra, using the virial method 
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Fig. 3. Comparison of the extended [O m] flux (again normalised 
to nuclear H/3) with two diff'erent measures of Fe ii strength. Left: 
Equivalent width of the Fe ii emission; right: Fe ii flux divided by 
the broad H/3 flux. Symbols as in Fig.|2] 



(IVestergaard & PetersonI 120061) . The method assumes that the 
broad line region (BLR) clouds are in gravitational equilibrium 
within the black hole potential, so that the central mass can be 
estimated as Mbh - fRv^/G where / is a geometry-dependent 
factor that needs to be calibrated empirically. The radius R of 
the BLR can be predicted from the quasar l uminosity using th e 
radius-luminosity scaling relation R oc Lp-^ JBentz et al.ll2006l) . 
The characteristic velocity of the emitting clouds can be esti- 
mated in two ways, from the FWH M of broad emission Unes or 
from the emission line dispersion cr. Collin et al] (12006') showed 
recently that the Mbh estimates are sensitive to the line profile 
shape and argued that criine gives more robust Mbh estimates 
than those based on FWHMiine. We compared both velocity mea- 
sures and found that while the differences were non-negligible, 
they are not important for the purpose of the present paper. We 
therefore a dopted the HP Line FWHM tog ether with the scaling 
relation by IVestergaard & PetersonI (|2006|) to obtain black hole 
masses. 

Having obtained black hole masses, we could also estimate 
the 'Eddington ratios' of each quasar, i.e. the actual accretion 
rates divided by the Eddington rates, or in terms of observables, 
the bolometric luminosities divided by the luminosities that each 
object would have if radiating at the Eddington limi t, Lboi/LEdd- 
We adopted a canonic bolometric correction (e.g., iKaspi et alj 
120001) of Lboi = 9/lL^(5100A). 



4. Discussion 

4.1. Relations between EELR detection and nuclear spectral 
properties 

In essentially all quasars where we detected extended emission 
in [Om], we also see a significant nuclear [Oiii] component. We 
also measured Fe ii emission to be very weak. This is shown in 
the left panel of Fig. |2] where we plot the nuclear [O iii]/Hy6 flux 
ratio against the Fe n equivalent width. Apparently, quasars with 
and without an EELR populate two distinct regions on this plot. 
By drawing a vertical line at log(Fe ii EW) = 1 .4, the two groups 
are well separated except for a single interloper with undetected 
EELR, which might be due to relatively poor observing condi- 
tion for that objects (seeing ~1.8")- Figure |3] compares two dif- 
ferent measures of relative Fe ii strength with the H/3 normalised 
flux ratio. In the left panel we employ the Fen EW as used al- 
ready in Fig. |2] and in the right panel we show the Fen/HyS 
ratio. While the objects with and without detected EELR are 
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Fig. 4. The EELR flux normalised by the nuclear H/3 flux plot- 
ted against continuum luminosity at 5 100 A (left), and against 
Eddington ratios (right). Symbols as in Fig.|2] 

again clearly separated in the former diagram, the distribution 
of Fe n/Hy6 values is essentially identical for the two groups. 

Comparing the normalised EELR flux with the FWHM of 
the broad H/? emission line we find an even clearer separation, 
albeit with a slightly larger number of interlopers. This is shown 
in the right-hand panel of Fig. |2] where 6/8 of all objects with a 
detected EELR have line widths well in excess of FWHM(Hy6) ~ 
4000 km/s. On the other hand, all quasars (except one) with non- 
detected EELR have 1000 km/s ^ FWHM(Hj8) ^ 3000 km/s. 
We note that the one exception, the RQQ PG 2214-1-139 with 
FWHM(H/?) ~ 5000 km/s and no detected EELR, shows a very 
peculiar multicomponent H/3 line shape. The two interlopers 
with FWHM(HyS) < 4000 km/s have the largest Fen EWs val- 
ues among the objects with detected EELR. 

4.2. Luminosities, black tiole masses and Eddington ratios 

We now explore whether the above results can be related to more 
fundamental quantities such as luminosities, black hole masses 
and Eddington ratios. Figure |4] presents, in a similar way as in 
Fig. |3] the EELR flux normalised by the nuclear H/3 flux plot- 
ted against the continuum luminosity (left panel) and against the 
Eddington ratio (right panel). All objects have broadly similar 
continuum luminosities, which basically ndes out the possibil- 
ity that it is simply the luminosity of the quasar determining 
the EELR properties. The FWHM is then the dominant term in 
the computation of Mbh, so that a pattern very similar to that 
seen in FWHM(Hj6) arises also for Mbh- We find that a dividing 
line at log(MBH/M0) = 8.5 separates the objects having smaller 
black hole masses and non-detected EELRs from the objects 
with larger masses and detected EELRs. If on the other hand we 
consider the Eddington ratios, there is again no clear separation 
between the two groups, with the distribution of log(Lboi/iEdd) 
covering the same range for detected and non-detected QSOs. 
We reiterate that this fact is independent of our choice of how to 
estimate black hole masses (Sect. [34] i. 

5. Conclusions 

By means of integral field spectroscopy, we have revealed the 
presence of extended emission line regions (EELRs) in a sam- 
ple of predominantly radio-quiet low-redshift QSOs. The sizes 
and luminosities are much higher than those of the well-known 
EELRs in Seyfert galaxies, extending out to several kpc and 
therefore stretching over the entire host galaxy (and possibly be- 
yond). Some 40 % of the objects in our sample show a prominent 
EELR, while we failed to detect any EELR signature in the re- 



maining 60 %. The detection or non-detection of an EELR can 
be linked to the nuclear spectral properties of the quasar itself, 
independen t of their radio c lassification, in excellent agr e emen t 
with what iBoroson & Ok3 (11984 and iBoroson et akl (1 19851) 
found from their ofi'-nuclear long-slit spectra. Our observations 
go further in that our integral field spectral data map out the en- 
tire nuclear-subtracted EELRs and derive their luminosities. 

Correlations between various nuclear s pectral properties 
were already studied by many authors (e.g. IBoroson & Greed 
119921; ISulentic et ai]|2000l; iNetzer et al.ll20()4l) . Here we provide 
new evidence that properties of the nuclear spectra, depending 
on the sub-pc scale conditions in AGN, are linked to properties 
of the host galaxies extending over several kpc. It thus appears 
possible to predict with high confidence whether or not a QSO 
will show a very extended emission line region just from its nu- 
clear spectrum: If it has broad H/? and weak Fe ii, it will proba- 
bly have not only strong nuclear, but also strong extended [O in] 
emission. 

It is tempting to interpret this relation between sub-pc and 
super-kpc scales in the context of the current discussion about 
AGN feedback affecting the host galaxies. Since we can rule out 
the QSO luminosity as a main driver of the differences between 
presence or absence of an EELR, the remaining possibilities are 
that there could be differences in the structure of the nuclear 
region (obscuration in connection with possible inclination ef- 
fects); or differences due to an outflow or jet, or the host galax- 
ies of the two QSO classes are intrinsically different, in particular 
with respect to their overall distribution of warm ionised gas. We 
conclude that QSO host galaxies come in at least two flavours: 
with or without extended ionized gas. Whether this highlights 
different initial conditions or a sequence of evolutionary stages 
we cannot say at present. More and better quality data will be 
needed to elucidate these intriguing trends. 
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